Gamma-Ray Constraint on Galactic Positron Production by MeV Dark Matter 
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The Galactic positrons, as observed by their annihilation gamma-ray line at 0.511 MeV, are diffi- 
cult to account for with astrophysical sources. It has been proposed that they are produced instead 
by dark matter annihilation or decay in the inner Galactic halo. To avoid other constraints, these 
processes are required to occur "invisibly," such that the eventual positron annihilation is the only 
detectable signal. However, electromagnetic radiative corrections to these processes inevitably pro- 
duce real gamma rays ( "internal bremsstrahlung" ) ; this emission violates COMPTEL and EGRET 
constraints unless the dark matter mass is less than about 20 MeV. 
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Introduction. — The SPI camera on the INTEGRAL 
satellite has recently observed the 0.511 MeV gamma- 
ray emission line arising from positron annihilation in the 
Galaxy [J. 0. Q- The flux from the Galactic center region 
is = (9.9^ x) x 10 -4 photons cm -2 s _1 , confirming 
earlier measurements H II @. The INTEGRAL data 
also provide new clues on the morphology of the emission 
region, which is consistent with a 2-dimensional gaussian 
of full width at half maximum (FWHM) of 9 degrees, 
with a 2 — (7 uncertainty range covering 6 — 18 degrees. 

A variety of astrophysical sources of positrons have 
been proposed, among them compact objects, massive 
stars, supernovae, gamma- ray bursts, and cosmic rays 
(see, e.g., Refs. pi m 111 USi llll)- However, these astro- 
physical sources have difficulty accounting for the inten- 
sity of the positron annihilation flux, and especially for 
the morphology of the emission region, which so far in 
the INTEGRAL data shows neither a disk component 
nor discrete sources. These facts motivate consideration 
of an exotic mechanism for the positron production, and 
one associated with the dark matter concentration at the 
Galactic center is naturally suggested. 

However, for dark matter candidates in the usually 
considered range of masses, 10 GeV-10 TeV, the produc- 
tion of positrons by dark matter annihilation would be 
accompanied by the production of other kinematically al- 
lowed particles (for a recent review of dark matter candi- 
dates, see Ref. |l2j). Even if direct production of gamma 
rays were suppressed, there would be a gamma-ray flux 
arising from the decays of secondaries (e.g., 7r°, which 
are created in the hadronization of quarks). For anni- 
hilation of typical dark matter candidates to account for 
the observed positrons, the associated flux of high-energy 
gamma rays in the direction of the Galactic center would 
exceed the EGRET data by several orders of magnitude. 

Boehm et al. p~3] have recently proposed that all of the 
characteristics of the observed signal could be well fit by 
a scenario in which light (1 — 100 MeV) dark matter parti- 



cles annihilate only into e + e~ pairs; the rate is controlled 
by the dark matter annihilation cross section (and its ve- 
locity dependence) , and the morphology by the assumed 
dark matter density profile pJJ LLJ, [l5j . Since the mass 
is below 100 MeV, the only other kinematically allowed 
annihilations would be to gamma rays and neutrinos. To 
avoid the direct gamma-ray constraints (among them the 
cosmic gamma-ray background pj|), and to account for 
the required positron production rates, these modes are 
postulated to not occur. While the assumed dark matter 
mass is quite low, it is claimed that this model is consis- 
tent with all laboratory data HI Q H [H [3 . In addi- 
tion, the model is apparently consistent with astrophysi- 
cal data and big-bang nucleosynthesis, depending on the 
assumptions made about the masses and couplings [lflj . 

In the model of Boehm et al., the positrons have a 
seemingly "invisible" birth, in that the dark matter anni- 
hilation only produces e + e~ pairs. Another clever aspect 
of this model is that positrons below 100 MeV will lose 
energy dominantly by ionization (avoiding production of 
possibly detectable inverse Compton and synchrotron ra- 
diation |2(j) and will remain confined to the Galactic cen- 
ter region (required to reproduce the observed emission 
region). The positrons thus remain invisible until their 
annihilation into 0.511 MeV gamma rays. Very simi- 
lar considerations apply to models of dark matter decay 
which produce the observed positrons 0, H^] • 

However, the dark matter annihilation process XX ~* 
e + e~ is necessarily accompanied by the process \X ~~ * 
e + e~7, arising from electromagnetic radiative correc- 
tions. This real gamma-ray emission, illustrated in Fig. 1, 
is known as internal bremsstrahlung, the name indicating 
that it arises from the Feynman diagram itself and is not 
due to propagation in a medium. In addition, the flux 
and spectrum of the internal bremsstrahlung gamma rays 
can be calculated with adequate accuracy without know- 
ing the new particle physics which mediates the dark 
matter annihilation. The internal bremsstrahlung gamma 
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FIG. 1: Diagram (a) shows the dark matter annihilation 
channel XX ~~* e + e~ assumed by Boehm et al. [l3j|: the 
positron production and propagation occur invisibly, with 
only the eventual positron annihilation being detectable. The 
open circle represents unspecified new physics. Diagram (b) 
shows the observation of this paper, that positron production 
must be accompanied by detectable internal bremsstrahlung 
(a similar diagram with radiation from the electron is not 
shown) . 



rays reveal the otherwise invisible dark matter annihila- 
tions, directly testing the central assumption of the pro- 
posed models, which is that the observed positrons were 
originally produced with energies up to 100 MeV. Astro- 
physical models would have to produce the same number 
of positrons, but at lower energies, typically a few MeV. 

We will show that for most of the proposed dark matter 
mass range, the flux of internal bremsstrahlung gamma 
rays would be inconsistent with COMPTEL and EGRET 
measurements of diffuse radiation from the Galactic cen- 
ter region, thus requiring the dark matter mass to be 
less than about 20 MeV. This result is almost completely 
independent of assumptions about the physical condi- 
tions at the Galactic center or the new particle physics 
which mediates the dark matter annihilation. In fact, 
it is more general than just dark matter annihilation or 
decay, and prohibits any model for the positron produc- 
tion in which the positrons are created at energies above 
about 20 MeV. We will also show that future Galactic 
gamma-ray data will be able to significantly improve the 
sensitivity of our constraint. 

Internal Bremsstrahlung. — The tree-level annihi- 
lation process XX ~~ * e + e~ will be subject to radiative 
corrections, which will affect the total cross section at 
0(a), where a — 1/137 is the fine-structure constant. 
These corrections are expected to be model-dependent 
but small. The portion of the radiative corrections that 
governs the emission of real gamma rays, i.e., the finite 
part of the outer radiative corrections, is also expected to 
be small, but is process-independent and factorizes from 
tree-level cross section (2^. Generally, those results ap- 
ply to soft gamma-ray radiation, but we will show that 
they are adequately accurate for our purposes, even for 
large gamma-ray energies. Internal bremsstrahlung has 
been observed in a very wide variety of processes (one 
with similar positron energies is muon decay HjJ). 

The process XX ~ ¥ e + e~ from x nearly at rest (typical 



halo velocities are 10 _3 c) produces monoenergetic elec- 
trons and positrons, with energies equal to the dark mat- 
ter mass m x . But when the gamma-ray energy is compa- 
rably large, it shares the phase space available to the final 
state. In order to reduce the process dependence, we use 
the final-state internal bremsstrahlung probability from 
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the kinematically similar process e + e 
modified by replacing the muon mass with the electron 
mass (since the initial state is XX-, it is appropriate 
to make this substitution while ignoring the identical- 
particle issue in e + e~ — > e + e~). 

The internal bremsstrahlung cross section is then 
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where E is the gamma-ray energy, and Otot is the tree- 
level cross section, which factors out, an essential point. 



Here s 

E, m e , 



4mr. and s' = Am x (m x — E). When m x ^> 



dE 



Otot x 



4a In (2m x /m e ) 
E 



(2) 



which displays the familiar scaling factors. The internal 
bremsstrahlung probability is very small for the gamma- 
ray energy range which we consider, which is well above 
the soft singularity. We find the most stringent con- 
straints when E is large, near m Xl so we need to account 
for the phase space corrections included in Eq. Q but 
not the simpler Eq. J2J). 

The flux of internal bremsstrahlung gamma rays is 
proportional to the positron production rate, which is 
determined from the 0.511 MeV intensity. Positron 
annihilation occurs either via direct annihilation into 
two 0.511 MeV gamma rays, or via the formation of a 
bound state called positronium. A singlet state (para- 
positronium) , which decays to two 0.511 MeV gamma 
rays, is formed 25% of the time, while the triplet state 
(ortho-positronium), which decays to three continuum 
gamma rays, is formed 75% of the time. The ratio of 
direct versus positronium annihilation can be measured 
by comparing the 0.511 MeV line intensity to the contin- 
uum intensity. It is customary to define the positronium 
fraction 1271 as 
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1.5 + 2.25($ 5U /$ C 



(3) 



This fraction depends on the physical conditions of the 
gas at the Galactic center, and observations suggest that 
/ = 0.93 ± 0.04, implying that most positrons annihilate 
via positronium |28j . This means that the dark matter 
annihilation rate and hence the internal bremsstrahlung 
rate are about 4 times larger than the positron annihila- 
tion rate observed via the 0.511 MeV line. 
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The spectrum of internal bremsstrahlung gamma rays 
per 0.511 MeV gamma ray is therefore 
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where we have accounted for the fraction of positrons 
that annihilate into 0.511 MeV gamma rays, and that 
each annihilation produces two 0.511 MeV gamma rays. 

Three Key Assumptions. — We now identify three 
key assumptions of the model of Boehm et al. [13] ; vio- 
lations of these assumptions act in the sense of strength- 
ening the internal bremsstrahlung constraint. 

First, it is assumed that the positron diffusion length is 
small compared to the size of the Galactic center region 
(this may not be the case; see Ref. ^3). Then the spatial 
distributions of dark matter annihilation and the internal 
bremsstrahlung gamma rays would be the same as the 
observed spatial distribution of e + e~ annihilation. If this 
assumption were violated, then the flux per steradian 
of internal bremsstrahlung gamma rays would be larger 
(and would come from a smaller angular region). 

Second, it is assumed that the positrons are brought 
to rest by ionization losses relatively quickly, so that the 
rates of dark matter annihilation and positron annihila- 
tion are in equilibrium. If this assumption were violated, 
then the positron annihilation rate today would reflect 
the dark matter annihilation rate in the past, while the 
internal bremsstrahlung rate would reflect the (larger) 
dark matter annihilation rate now. 

Third, it is assumed that the positrons annihilate at 
rest. If this assumption were violated, then to pro- 
duce the same 0.511 MeV flux, the dark matter annihila- 
tion and internal bremsstrahlung gamma-ray rates would 
have to be larger. Additionally, positron annihilation in 
flight could produce a detectable high-energy gamma-ray 
signature of their own, though it would depend on the 
details of the physical conditions in the Galactic center. 

COMPTEL/EGRET Constraints. The flux of 
internal bremsstrahlung gamma rays per steradian is 
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where we have assumed that 1/2 of the 0.511 MeV flux is 
emitted from an angular region AQ given by the gaussian 
FWHM of 9 degrees. This becomes 
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The corresponding internal bremsstrahlung gamma-ray 
spectra are shown in Fig. 2. The energy dependence of 
the cross section in Eq. (JJ) gives spectra which peak 
at E ~ m x , and fall off as 1/E; for easier compari- 
son with the data, we plot E 2 times the spectra. The 
offsets between the solid lines at low E reflect just 
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FIG. 2: The solid lines show the internal bremsstrahlung 
spectra labeled by the assumed dark matter masses in MeV, 
and normalized to produce the observed Galactic 0.511 MeV 
flux (with the light dotted lines corresponding to its uncer- 
tainty); we assumed that the FWHM of the 0.511 MeV emis- 
sion region is 9 degrees. The approximate COMPTEL (a) 
and EGRET (b) data and their uncertainty band are shown 
by the box, and the long-dashed line below it indicates the 
rough constraint on the maximum allowed contribution of in- 
ternal bremsstrahlung gamma-rays. 



the ln(2m x /m e ) factor in the approximate and process- 
independent Eq. (J2J. The mild process dependence (the 
difference between Eq. Q and Eq. in the treatment 
of the phase space) is seen in Fig. 2 in the slight turnovers 
in the solid lines before the endpoints, as well as a gen- 
eral reduction by a factor < 2. Had we used just Eq. J5J, 
our constraints would have been somewhat too strong. 

We can constrain the internal bremsstrahlung contri- 
bution by comparing to COMPTEL and EGRET mea- 
surements of the diffuse gamma-ray flux from the Galac- 
tic center region |29t yCJ, |31| . Those fluxes were averaged 
over a Galactic center region of ±5 degrees in latitude and 
±30 degrees in longitude, but were shown to have mild 
variation across that region. The observed 0.511 MeV 
emission is from a smaller region of FWHM 9 degrees. 
Thus in order to constrain an internal bremsstrahlung 
contribution to the central circle, we must first subtract 
the diffuse astrophysical contribution as measured over 
the whole COMPTEL/EGRET rectangle. We have con- 
servatively assumed that the COMPTEL/EGRET data 
band is given at 1 sigma, and that an excess contribution 
of more than about 50% could not be tolerated. A fit to 
the full energy and angular dependence of the signal and 
background would yield a more stringent constraint. 

Even with our very conservative treatment, we can eas- 
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ily constrain the dark matter mass in the Boehm et al. 
mo del Jl3l| to m x < 20 MeV (for decaying dark matter 
|2lL |22| we determine a similar upper limit). The un- 
certainties on the COMPTEL/EGRET data (as well as 
our approximate handling thereof), the 0.511 MeV flux, 
and the positronium fraction are relatively unimportant, 
though improvements would be welcomed. The largest 
uncertainty on the constraint arises from the uncertainty 
in the size of the 0.511 MeV emission region. Following 
Refs. 0, 0, EJ El, we assumed 9 degrees, but the 2 — a 
range spans 6 — 18 degrees, allowing an order of magni- 
tude in the internal bremsstrahlung flux per steradian; 
this corresponds to varying the upper limit on m x from 
10 to 60 MeV. Improved results from INTEGRAL on the 
size of the emission region are thus eagerly awaited. 

Discussion and Conclusions. — In order to explain 
the Galactic positron excess, as well as the smooth and 
centrally symmetric morphology of the 0.511 MeV emis- 
sion observed by. INTEGRAL^, H , Boehm et al. 
(see also Refs. [lj,llall3ll3ill3) proposed an intriguing 
model where e + e~ pairs are produced by the annihilation 
of light (1 — 100 MeV) dark matter candidates, which 
apparently can evade all present accelerator and astro- 
physical constraints. We have pointed out the central as- 
sumption of this model, i.e., that energetic positrons may 
be produced in the Galactic center with no other obser- 
vational consequences other than their eventual annihi- 
lation into 0.511 MeV gamma rays. We have shown that 
such an "invisible" birth is prohibited by the emission 
of internal bremsstrahlung gamma rays from the original 
dark matter annihilations unless the dark matter mass 
is less than about 20 MeV, disallowing most of the pro- 
posed range. We have arrived at this constraint in a very 
conservative fashion, and expect that improved gamma- 
ray data and a more sophisticated analysis will signif- 
icantly improve the sensitivity. Our constraint is very 
nearly independent of the new physics of the dark mat- 
ter particles, and would be strengthened by relaxing the 
assumptions on the physical conditions at the Galactic 
center. 

Although our results have been presented in the con- 
text of dark matter annihilation, they are much more gen- 
eral. Any mechanism which produces energetic positrons 
will be accompanied by internal bremsstrahlung gamma 
rays, and the rate is nearly independent of the tree level 
cross section. Thus for any mechanism that creates 
enough positrons to account for the 0.511 MeV line, if 
those positrons are produced above about 20 MeV, the 
accompanying internal bremsstrahlung will violate the 
COMPTEL/EGRET constraints. Dark matter decay is 
an example of such a mechanism [^, ■ The ultimate 
sensitivity of the technique at low energies will be reached 
when the dominant sources are astrophysical sites which 
produce gamma rays and positrons at comparable (tree- 
level) rates. 
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